High-field phase diagram of the Haldane-gap antiferromagnet Ni(C5Hi4N2)2N3(PF6) 
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We have determined the magnetic phase diagram of the quasi-one-dimensional S = 1 Heisenberg 
antiferromagnet Ni(C5Hi4N2)2N3(PF6) by specific heat measurements to 150 mK in temperature 
and 32 T in the magnetic field. When the field is applied along the spin-chain direction, a new phase 
appears at Hc2 ~ 14 T. For the previously known phases of field-induced order, an accurate deter- 
mination is made of the power-law exponents of the ordering temperature near the zero-temperature 
critical field He, owing to the four- fold improvement of the minimum temperature over the previous 
work. The results are compared with the predictions based on the Bose-Einstein condensation of 
triplet excitations. Substituting deuterium for hydrogen is found to slightly reduce the interchain 
exchange. 

PACS numbers: 75.30.Kz, 75.40.Cx, 75.50.Ee 



I. INTRODUCTION 

One-dimensional (ID) integer-spin Heisenberg antifer- 
romagnets (HAF) are well-known for the Haldane energy 
gapi between the singlet spin-liquid ground state and the 
lowest excited state, which is an S* = 1 triplet. The ap- 
plication of a magnetic field leads to Zeeman splitting 
of the triplet and eventual vanishing of the gap A at 
He ^ A/g/iB, where the energy of the lowest branch of 
the split triplet reaches the ground-state level. At this 
critical field, the expected quantum transition is to the 
Tomonaga-Luttinger spin liquid)^ in which spin correla- 
tions decay with characteristic power laws. This scenario 
remains robust in an array of ID integer-spin HAF chains 
against the introduction of interchain exchange. Such a 
coupling reduces He but, provided it is small, does not 
destroy the singlet ground state below He- Above He, 
it leads to finite-temperature long-range order (LRO), 
which can be described as the three-dimensional Bose- 
Einstein condensation of the lower-branch tripletsA^iSiLS 
The nature of the ordered state, including the robustness 
of the Tomonaga-Luttinger spin liquid, is of strong cur- 
rent interest. 

The first experimental evidence^ of a Haldane gap 
was found in CsNiCls. This material has a relatively 
large interchain couphng: J'/ J — 0.017, where J' 
and J are interchain and in-chain exchanges, respec- 
tively. As a consequence, Neel ordering occurs at 4.85 K. 
Ni(C2H8N2)2N02(C104) (NENP), with the smaU \ J'\/J 
ratio of 4 X 10"*, was the first Haldane-gap antiferro- 
magnet with no Neel orderingiiSiii Magnetic suscepti- 
bility shows no anomaly indicative of ordering, at least 
down to 4 mK.^^ Naturally, this material was a promis- 
ing candidate for the field-induced LRO at fields above 
He- However, it has been found that such an order is 
preempted by the staggered effective field arising from 
the staggered g tensors of the Ni^+ ions within each spin 
chain. 



To date, another spin-1 chain material 
Ni(C5Hi4N2)2N3(PF6) (NDMAP) remains the only 
laboratory model of a ID HAF array in which the nature 
of the field-induced LRO states has been revealed by 
experiment. ^'^ This distinction owes to the unusually low 
He, which is readily accessible to many experimental 
probes, as well as the absence of a staggered field. 
NDMAP has an orthorhombic structure with the lattice 
parameters a = 18.046 A, b = 8.7050 A, and c = 6.139 A, 
with the antiferromagnetic spin chains running along 
the c axis»i^ According to inelastic neutron scatteringfi^ 
the in-chain exchange is J = 33.1 K, and the relative 
strengths of the interchain exchanges are J^/ J « 6 x 10"'' 
and \Ja\/J < 10"'' along the b and a axes, respectively. 
Easy- plane crystal-field anisotropy where z is the 

crystallographic c axis and D/J = 0.25^^ is responsible 
for the strong anisotropy of both the gap A and the 
magnetic phase diagram of the field- induced LRO states, 
as has been measured by magnetic susceptibility^ 
specific heatfiiiSiii magnetization,iSiia ESR^SfliSi and 
neutron scatteringi^^i^'^i^^-^'^'^"' 

The gap energies of the triplet, whose degeneracy is 
lifted by the crystal-field anisotropy, have been measured 
by inelastic neutron scattering^^ to be Aq = 0.42 meV, 
Afc = 0.52 mcV, and A^ = 1.9 meV for the Si = 
excitations quantized along the a, b, and c axes, respec- 
tively. Here z is the quantization axis. The small split- 
ting of a doublet into A^ and A^ is due to a weak in-plane 
anisotropy E{S'^ — Sy). Because of the gap anisotropy. He 
depends on the crystal orientation and ranges from about 
4 T to 6 TiSiiSJi Above He, the ordered states are com- 
mensurate with the crystal lattice for all four magnetic- 
field directions studied by neutron diffractioniiiSi For 
field H applied along the a axis and along [110], the or- 
der is short-ranged and confined within each be plane, 
whereas it is long-ranged and three-dimensional for H\\c 
and ir||[011]. 

The neutron experiments have detected no incommen- 
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surate modulation of the magnetization component par- 
allel to the field. Such a modulation has been predicted 
for the Tomonaga-Luttinger spin liquid^S*^ and its ab- 
sence in NDMAP can be understood as a consequence 
of a lack of axial symmetry, which is necessary for the 
existence of the Tomonaga-Luttinger spin liquidi^ It is 
important to note that even H\\c is not an axisymmetric 
field configuration for NDMAP because of its unique ge- 
ometry. In this material, the crystal-field anisotropy of 
Ni^+ is determined by the local symmetry of the NiNg 
octahedron. First, there is a weak in-plane anisotropy 
as described above. Second and probably more impor- 
tant, the principal axis of the octahedron is tilted from 
the c axis by 15° toward the a axis, with the tilt direc- 
tion alternating from chain to chaiufM Therefore, there 
is no magnetic-field direction that strictly satisfies axial 
symmetry. 

In this paper, we explore the field-induced LRO states 
to 32 T in field and 150 mK in temperature by means of 
specific heat measurements. The present study greatly 
extends the H-T parameter space for the phase diagram, 
which has been previously limitedi^'iS'ii to 12 T in H 
and 0.52 K in temperature T except for a few isolated 
points obtained by neutron diffractioni2*24 at lower tem- 
peratures with less precision. Preliminary results have 
been reported in Ref. |2^ 



o 
E 



O 




T(K) 



0.4 



0.3 



0.2 



0.1 



0.0 



(b) 



T ^ 

1—0.83 K 
— o— 0.77 K 
— i— 0.71 K 
—a— 0.66 K 
—♦—0.59 K 
— o— 0.51 K 
— T— 0.45 K 
— V— 0.35 K 
—■—0.30 K 
—□—0.22 K 
— ®— 0.15 K 



/ h ° 



H//C 



4 5 
H(T) 



II. EXPERIMENTAL PROCEDURES AND 
RESULTS 

The single crystals of NDMAP used in this work were 
grown from aqueous solutions by the method described in 
Ref. 0. Fully deuterated crystals were used to eliminate 
the nuclear specific heat of protons in the set of measure- 
ments extending to the lowest temperature, whereas a 
hydrogenous sample was used when adequate. The spe- 
cific heat measurements to 18 T were performed in a su- 
perconducting magnet at temperatures down to 150 mK 
with a dilution refrigerator. The relaxation calorimeter 
for this setup has been described in Ref.lU For measure- 
ments from 20 T to 32 T, another relaxation calorimeter 
with a built-in ^He refrigerator was used in a resistive 
Bitter magnet. Each sample for the H\\c configuration 
was attached with silver pain^ to the vertical face of 
a small silver bracket, whose horizontal face was in turn 
glued with a Wakefield compound^^ to the calorimeter 
platform. The sample for the H\\a configuration was 
mounted with silver paint on a piece of 0.13 mm-thick 
sapphire substrate, which in turn was glued with a Wake- 
field compound to the platform. A total of three samples 
ranging in mass from 4 mg to 8 mg were used to cover 
the different field ranges and orientations. 

The magnetic- field and temperature dependence of the 
specific heat of a deuterated sample is shown in Fig. 1 
for fields up to 18 T, applied along the c axis of the crys- 
tal. At each field, the peak in the specific heat seen in 
Fig. 1(a) clearly indicates a phase transition. As can be 



FIG. 1: Specific heat of deuterated NDMAP for a magnetic 
field parallel to the c axis, (a) Temperature dependence at 
constant fields. The inset gives an expanded view of the region 
near the specific heat peak, (b) Magnetic-field dependence at 
constant temperatures. The lines are guides for the eye. 



seen in the inset, the transition temperature denoted by 
the peak position first increases with increasing field up 
to about 10 T and then decreases for fields up to 14 T, 
where it starts to increase again, making a shallow local 
minimum at around 14 T. Furthermore, the peak height 
has a less pronounced minimum at roughly the same field. 
These features reveal the existence of a new phase bound- 
ary, which separates two field-induced ordered states at 
around 14 T. For temperatures less than 0.85 K, spe- 
cific heat is shown as a function of the magnetic field 
in Fig. 1(b). Again, the peak clearly indicates a phase 
transition. We have confirmed that the transition tem- 
peratures determined by such magnetic-field scans agree 
with those determined by the temperature scans. 

The specific heat for the same field direction at higher 
fields produced by the resistive magnet are shown in 
Fig. 2. These data were obtained with the hydrogenous 
sample. At 26 T, the nuclear contribution of the protons 
is visible at temperatures below 0.7 K. However, the peak 
due to transition stands out, since it occurs at a higher 
temperature. In this field region extending from 20 T to 
32 T, the transition temperature obtained from the peak 
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FIG. 2; Specific lieat of liydrogenous NDMAP at constant 
magnetic fields ranging from 20 T to 32 T, applied along tlie 
c axis. Tlie lines are guides for the eye. 
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FIG. 4: Magnetic-field dependence of the specific heat of 
deuterated NDMAP measured at constant temperatures. The 
magnetic field was applied along the a axis. The lines are 
guides for the eye. 
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FIG. 3: Magnetic-field dependence of the specific heat of 
hydrogenous NDMAP in the ordered phases. The field was 
applied along the c axis. The anomaly at the lowest field 
of each curve is due to the proximity of the transition to 
the disordered spin-liquid phase. The uninteresting proton 
contribution raises the curve at 0.26-0.29 K with respect to 
those at higher temperatures. The lines are guides for the 
eye. 



position varies only monotonically, and so does the peak 
height. 

To further investigate the phase diagram near 14 T, we 
have measured the specific heat of the same hydrogenous 
sample as a function of the magnetic field, again applied 
along the c axis, but to higher fields than in Fig. 1(b). In 
these measurements, a constant electric current was fed 
to the heater of the thermal reservoir of the calorimeter, 
allowing the temperature to rise monotonically with an 
increasing field as dictated by the magnetoresistance of 
the heater. As seen in Fig. 3, a plateau-like anomaly in 



the specific heat occurs at around 14 T, clearly indicat- 
ing a phase transition. This anomaly was overlooked in 
our preliminary report, where the field range investi- 
gated was too narrow. Within the experimental resolu- 
tion, no corresponding feature is found in magnetization, 
which has been measured at 80 niK and 1.3 K using a 
pulsed magnet and is featureless up to 60 T except for 
an anomaly at 

For magnetic fields applied along the crystallographic 
a axis, which is perpendicular to the spin chains, another 
deuterated sample was used to measure the specific heat 
at temperatures below 2.7 K as a function of field, as 
shown in Fig. 4. Again, a transition to the field-induced 
ordered phase is clearly indicated by a peak at each tem- 
perature. For fields ranging from 6.2 T to 6.5 T, these 
measurements were supplemented by temperature sweeps 
similar to those shown in Fig. 1(a). No new phase bound- 
ary was found for this field direction up to 18 T. 

The magnetic phase diagram determined from the po- 
sitions of the specific-heat peaks and the plateau-like 
anomalies is shown in Fig. 5 for the two magnetic-field 
directions studied. When the field is along the c axis, the 
transition temperature exhibits a shallow but distinct lo- 
cal minimum at about Hc2 ~ 14 T, and a new phase 
boundary extends nearly horizontally from this mini- 
mum. The small break in the phase boundary between 
the new high-field phase and the thermally disordered 
phase at around 20 T indicates that the hydrogenous 
sample has a slightly higher transition temperature than 
the deuterated sample. As stated earlier, the calorime- 
ters used for the two samples were different. However, we 
have confirmed in another experiment the consistency of 
the temperature scales of the two calorimeters. There- 
fore, the discontinuity in the phase boundary is not an 
artifact and indicates that deuterated NDMAP in fact 
has a somewhat weaker exchange. When the field is along 
the a axis, the transition temperature rises rapidly with 
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FIG. 5: Field-temperature phase diagram of NDMAP. The 
open symbols are for the field applied along the c axis and 
the closed symbols for the field along the a axis. The circles 
are for the deuterated samples, and the squares and diamonds 
for the hydrogenous sample. The lines are guides for the eye. 



FIG. 6: Field-temperature phase diagram of NDMAP below 
20 T for the field applied parallel to the c and a axes. The 
closed symbols are for the deuterated samples of the present 
work, where the circles are from the temperature-scan data 
and the diamonds from the field-scan data. The open circles 
are for hydrogenous samples from Ref. Il3l The lines are low- 
temperature fits of the phase boundaries of the deuterated 
samples to the expression Tc (x (H — He)" ■ 



increasing field, with no indication of a new phase. 

III. DISCUSSION 

This is the first time more than one field-induced or- 
dered phase has been found in a quasi-lD antifcrromag- 
net for a given magnetic-field orientation. We believe 
it is significant that the field direction required for the 
new phase is along the crystalline c axis. This direc- 
tion, parallel to the spin chains, is quite unlike the a, 
b, and [Oil] directions in that the ordering field rises 
very rapidly with temperature in the overall H-T phase 
diagramiiiiLii This feature probably exemplifies the un- 
derlying tendency of the spin chains in this field configu- 
ration to form a Tomonaga-Luttinger spin liquid, as has 
been presumed earlier. Hence it is likely that identifying 
the new phase above Hc2 and understanding the mech- 
anism of the transition will shed light on the physics of 
the Tomonaga-Luttinger spin liquid. 

At present, however, with the absence of experimen- 
tal data in the field region near and above Hc2 other 
than the specific heat and magnetization, which shows 
no anomaly at Hc2, we can at best speculate on the 
nature of the transition. One possible scenario can be 
a transition involving spin rotation around the direc- 
tion of the applied field, which is parallel to the spin 
chains. Such an exotic transition has been observed for 
instance in a spin-1/2 chain material BaCu2Si207, where 
it is presumably driven by a competition between an off- 
diagonal exchange and the magnetic fieldi^ Although 
the Dzyaloshinskii-Moriya interaction^i the likely candi- 
date for the off-diagonal exchange, is usually quite small 
for Ni^+ and has not been detected in NDMAP, a spin- 



rotation transition is an interesting possibility. 

Recently, WangS^ has considered a spin-1 chain with 
broken symmetry using a fermionic field theory and has 
predicted that a magnetic field larger than He will re- 
store an approximate axial symmetry and lead to a well- 
defined second transition from a commensurate phase 
to an incommensurate phase. Unlike in the Tomonaga- 
Luttinger spin liquid, the excitations in the incommen- 
surate phase will be gapped. However, the gap will be 
small and roughly equal to (Af, — Aq)/2, when the field is 
applied along the chain direction c. Although Wang has 
predicted an absence of anomaly in specific heat — as 
well as magnetization — at the second transition, which 
is other than first-order, it is worth exploring the possi- 
bility of the phase above Hc2 being an incommensurate 
phase. 

The Bose-Einstein condensation of triplets is believed 
to be a valid description of field-induced LRO in all 
gapped antiferromagnets, regardless of the type of dis- 
ordered ground state below He and the type of the or- 
dered state above it, provided that the lowest excitations 
are bosons. According to the Bose-Einstein condensa- 
tion theory^i^ the phase boundary in the H-T param- 
eter space obeys a power law Tc cx {H — He)"- The 
field-induced ordering in the 5 = 1/2 spin-dimer material 
TlCuCls was the first case that was analyzed in terms of 
the Bose-Einstein condensationiS However, the exponent 
a of 0.50 determinedSi for this material is smaller than 
the theoretical value of 2/3 obtained by a Hartree-Fock 
approximation for a quadratic dispersion for the triplets. 
The exponent for Ni(C5Hi4N2)2N3(C104) (NDMAZ), a 
spin-1 chain material'^'* similar to NDMAP, has been re- 
ported to be 0.45 for the magnetic field applied along the 
c axis.'^^ This value is also significantly smaller than the 
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theoretical one. 

The four-fold improvement in the minimum tempera- 
ture over the previous worl*iii& for the field directions 
along the c and a axes allows us to determine the field 
dependence of the ordering temperature near He with 
high accuracy. The relevant portions of the phase bound- 
aries from Fig. 5 are reproduced in Fig. 6, separately 
for the two field directions. These boundaries, shown 
with closed symbols, are for deuterated samples. Excel- 
lent fits of the data below 0.77 K for H\\c and below 
1.9 K for H\\a are obtained with the power law, giving 
a = 0.34 ± 0.02 and 0.57 ± 0.01, respectively. The crit- 
ical fields are = 3.82 ± 0.03 T and 6.01 ± 0.04 T, 
respectively, along the c and a axes. 

Recently, specific heat measurements of the 
S = 1 bond-alternating-chain antiferromagnet 
Ni(C9H24N4)N02(C104) (NTENP) have found 
a = 0.334 and 0.52, respectively, for the field par- 
allel and perpendicular to the spin chains. "^^ The 
similarity of these exponents to those for NDMAP 
suggests that NTENP has similar field-induced ordered 
phases, although its low-field ground state is a spin-dimer 
singlet'^^ rather than a Haldane spin liquid. 

Including ours, all the experimental values found to 
date for the exponent a are smaller than 2/3 pre- 
dicted by the theory. However, the values for NDMAP 
and NTENP for the magnetic field applied along the 
chain direction agree, within the combined uncertain- 
ties of the experiments and the calculation, with the 
quantum Monte-Carlo results by Wessel et al.r^ who 
have found a = 0.37 ± 0.03 for S = 1/2 spin dimers 
with a weak three-dimensional inter-dimer coupling and 
a — 0.32 ± 0.03 for two-leg spin-1/2 ladders with a 
weak three-dimensional inter-ladder coupling. Their re- 
sults support the Bose-Einstein condensation theory, ac- 
cording to the quantum Monte-Carlo calculation for 
three-dimensionally coupled 5 = 1/2 spin dimers by 
Nohadani et al.,^^ who have shown that a deviates 
downward from 2/3, as the temperature range of the 
power-law fit is widened. Furthermore, Misguich and 
Oshikawa^" have used a realistic dispersion for the 
triplets in TlCuCla in their calculation of the field de- 
pendence of the Bose-Einstein condensation temperature 
Tc and found good agreement with the experiment.® It 
remains to be seen, however, whether the Bose-Einstein 
condensation theory can explain the strongly anisotropic 
a for NDMAP and NTENP, particularly the extremely 
small a when the field is applied along the spin chains. 

In Fig. 6, we have included the phase boundaries for 
hydrogenous samples from the previous work, as shown 
with open symbolsiifiJ'^ Again, the ordering temperature 
is somewhat higher for the hydrogenous samples except 
at the fields of 10 T and 12 T applied parallel to the c axis, 
whereas He is identical for a given field direction within 
our accuracy. This indicates that the interchain exchange 
J' is slightly weaker in deuterated NDMAP than in hy- 



drogenous NDMAP, since decreasing J' lowers the order- 
ing temperature directly but has little effect on H^, which 
is primarily determined by the in-chain exchange J. A 
previous study, ''^ which compared the phase diagrams of 
deuterated and hydrogenous samples for H along the b 
axis, saw a much smaller but qualitatively similar differ- 
ence. 

We propose that the smaller J' for deuterated NDMAP 
is caused by the smaller zero-point motion of deuterons 
leading to less overlap of the electronic wavefunc- 
tions in the exchange paths involving hydrogens. In 
some classical antiferromagnets, such as the linear- 
chain compound (CIl3)4NMnCl3 (TMMC), deuteration 
has been reported to reduce the Neel temperaturei^Siii 
In the quasi- two-dimensional organic salt k-(BEDT- 
TTF)2Cu[N(CN)2]Br, substituting deuterium for the hy- 
drogen of the ethylene groups has been shown to drive 
the system from a superconductor to an antiferromag- 
netic Mott insulator44 To our knowledge, however, the 
present work represents the first observation of an effect 
of deuteration in a Haldane-gap antiferromagnet. 

In summary, the magnetic phase diagram of NDMAP 
has been extended by specific-heat measurements to 
150 mK in temperature and 32 T in the magnetic field. A 
new transition has been found at Hc2 ~ 14 T, when the 
field is parallel to the spin chains. Further study using 
techniques other than specific heat is needed to investi- 
gate the nature of the new phase, including the possibil- 
ity that it is an incommensurate phase. In addition, we 
have determined the power-law exponents for the transi- 
tion temperatures of the field-induced ordering from the 
spin-liquid state and compared them with the exponents 
for TlCuCla, NDMAZ, and NTENP and with the predic- 
tions of the Bose-Einstein condensation theory. Finally, 
we have observed an effect of deuterium substitution on 
the ordering temperature. The result indicates that the 
substitution slightly decreases the interchain exchange J' 
but hardly affects the in-chain exchange J. 
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